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£>y Abstract 

a 

The author considers an elliptic analogue of the Knizhnik-Zamolodchikov equa- 
tions - the consistent system of linear differential equations arising from the elliptic 
solution of the classical Yang-Baxter equation for the Lie algebra si n ■ 

The solutions of this system are interpreted as traces of products of intertwining 
operators between certain representations of the afhne Lie algebra sIn- 

A new differential equation for such traces characterizing their behavior under 
the variation of the modulus of the underlying elliptic curve is deduced. This 
equation is consistent with the original system. 

It is shown that the system extended by the new equation is modular invari- 
' ant, and the corresponding monodromy representations of the modular group are 

^\ , defined. 

Some elementary examples in which the system can be solved explicitly (in terms 
of elliptic and modular functions) are considered. Another example leads to a spe- 
cial case of Heun's equation - a second order Fuchsian equation with four singular 
<D , points. 

The monodromy of the system is explicitly computed with the help of the trace 
interpretation of solutions. Projective representations of the braid group of the 
torus and representations of the double affine Hecke algebra are obtained. 
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5t i Introduction 

In 1984 Knizhnik and Zamolodchikov [KZ] studied matrix elements of products 
of intertwining operators between representations of the affinization g of a finite 
dimensional simple complex Lie algebra q at level k. These matrix elements are 
analytic functions of several complex variables, and it was found that they satisfy 
a certain remarkable system of linear differential equations which is now called the 
Knizhnik-Zamolodchikov (KZ) system: 

3 = 1,3^1 

Here ^(^i, ■■■,z n ) is a function of n complex variables with values in the product 
W = Vi (S> V2 <8> • • • <8> V n of n representations of g, h y is the dual Coxeter number 
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of g, and 
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(2) 



(x p )i — Idi <g> • • • <g> Idi_i <g> x p <g> Id l+ i ® • • • (g) Id„ e 



End(VK) 



where the summation is over an orthonormal base {x p } of g with respect to the 
invariant form, and Idj denotes the identity operator in Vj . 

Solutions of the KZ equations are very interesting special functions which gener- 
alize the Gauss hypergeometric function. By now an explicit integral represen- 
tations of these functions has been found [M],[SV], and their monodromy has 
been completely computed [Koh] , [Dr] , [TK] (for g = sZ 2 ),[V]. It turned out that 
the monodromy of the KZ system expresses in terms of the quantum i?-matrix 
- the quasitriangular structure of the corresponding quantum group U q (g), where 



Cherednik [Chi] considered a general consistent system of differential equations 
of the "factorized" form 



where r(u) is a meromorphic function with values in g ® g, and rij(u) denotes the 
action of r(u) in W: the first factor acts in Vi and the second one in Vj, and the 
unitarity property Tij{u) = —rji(-u) is assumed (this condition is equivalent to 
the invariance of solutions of system (3) under the simultaneous translation of all 
variables zj by the same constant). It was observed in [Chi] that system (3) is 
consistent if and only if the function r(u) is a classical r-matrix, i.e. if it satisfies 
the classical Yang-Baxter equation: 



[rij(zi - zj), rik{zj - Zk)] + [nj(zj - z k ),r jk (z k - Zj)] + [r lk (z k - z l ),r jk (z i - zj)] = 0. 

Therefore, a consistent system (3) is called a local r-matrix system. 

Clearly, the KZ system is a special case of a local r-matrix system, for a simplest 
r-matrix r(u) = 0,/u, O = J2 V 

(the summation is over an orthonormal base 
{x p } of g). This gives rise to a question: what other r-matrix systems are possible? 
This question was essentially answered by Belavin and Drinfeld in 1982 [BeDr]. 
They classified all solutions of (3) satisfying the nondegeneracy condition: r(u) is 
invertible as a map g* — > g for at least one complex number u. This classification 
states that all such solutions are unitary and, in terms of dependence on u, there 
are only three types of functions r(u): rational, trigonometric, and elliptic. For 
instance, the function r{u) = Cl/u which is involved in the KZ system is a rational 
r-matrix. 

Trigonometric an elliptic nondegenerate unitary solutions of the classical Yang- 
Baxter equation are completely classified [BeDr]. On the contrary, a satisfactory 
classification of rational solutions is unknown. 

The Belavin-Drinfcld classification suggests a two-step generalization of the KZ 
system: KZ equations (rational r-matrix equations) - trigonometric r-matrix equa- 
tions elliptic r-matrix equations. One should expect that these local r-matrix 
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systems should have remarkable properties and provide new interesting special func- 
tions as their solutions. 

Cherednik [Ch2] found an interpretation of solutions of nondegenerate unitary 
r-matrix equations in terms of representation theory of affinc Lie algebras. It was 
proved in [Ch2] that the general solution is the coinvariant (or the so-called r- 
function) of the Lie algebra of g-valued rational functions on a rational or elliptic 
curve with singularities at designated points Z\, z n in a certain representation of 
this algebra. This interpretation has found many applications. 

The properties of the trigonometric r-matrix equations are now fairly well un- 
derstood. Let g = n + f) n~ be a Cartan decomposition of g, and let Cl + , £l~ , 
Cl° be the orthogonal projections of O to the subspaces n + ® n~, n~ <£> n + , and 
t) (g> f), respectively . Then the simplest trigonometric g <g> g- valued r-matrix has 
the form r{u) = ° + „_ 4 L ° , where = + ^Qq. It turns out that there exists 
a transformation of coordinates which maps the n + 1-point KZ equations to the 
n-point trigonometric r-matrix equations with the above r(u) (which are called the 
trigonometric KZ equations), which allows one to carry over to this case all the 
results about the KZ equations, including the integral formulas. 

Trigonometric r-matrix equations with more general r-matrices were studied by 
Cherednik [Ch2], who gave an explicit integral formula for the general solution. 

Nondegenerate unitary elliptic r-matrices exist only for g = slj^. They were 
found in the case N — 2 by Sklyanin and in the general case by Belavin[Be] (note 
that the N = 2 elliptic r-matrix is the quasiclassical limit of Baxter's quantum 
i?-matrix which arises in statistical mechanics): 
(5) 

P (z\r) = n((z)+ ]T (i^« 7 -)(o)( C (,_!^| T )_ C (!^| T )), 

a<m,n<N-l,rn 2 +n 2 >0 

where /3, 7 are two commuting inner automorphisms of sIn of order N with no 
common invariant vectors, r lies in the upper half of the complex plane, and ((z\t) is 
the Weierstrass ^-function. The properties of the corresponding r-matrix equations 
are not very well understood. For example, integral formulas (or any other explicit 
representations) for solutions are unknown. 

It has been anticipated that some progress in the elliptic case can be achieved 
by using the intertwining (vertex) operator language which was originally used 
by Knizhnik and Zamolodchikov[KZ]. Frenkel and Reshetikhin [FR] conjectured 
that if one takes traces of products of intertwiners rather than matrix elements, 
one should be able to obtain solutions to the elliptic r-matrix equations. The 
same idea occurs in the paper of Bernard [Ber] who studied expectation values 
for the Wess-Zumino-Wittcn model on an elliptic curve and was led to consider 
trace expressions of a similar sort. Bernard deduced some differential relations for 
traces, but they were not a closed system of differential equations since the two 
commuting inner automorphisms did not enter the game. The idea to consider 
traces is also suggested by the aforementioned theory developed in [Ch2] which 
states that solutions of an elliptic r-matrix system should express in terms of the 
r-function for the corresponding elliptic curve, which is basically a trace expression 
involving vertex operators. 

This paper is devoted to making this idea work (at the price of a few modifica- 
tions). More precisely, we represent solutions of the elliptic r-matrix equations in 
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the form 

(6) F(z u z n \q) = Tr \ MKk (<P\ Zl ) . . . $ n (z n )Bq- 9 ), 

where is a Verma module over (a twisted version of) the affinc algebra sIn, 

&(zi) : M\ it k — > M\ i _ 1 ,k ® Vi are intertwiners for s^at (the action of sZjv on V*: 
(a<8)t m )u = z m av, a E g, v e V*; the hat denotes a completion of the Verma module), 
q G C, |g| < 1, 9 is the grading operator in M^fc, and B is the map of Verma 
modules induced by an outer automorphism of sIn of order N (rotation of the 
affine Dynkin diagram, which is a regular 7V-gon, through the angle of 2n/N). This 
representation helps to compute the monodromy of the elliptic r-matrix equations. 
It is still unclear how to deduce integral formulas for solutions similar to those 
existing for the rational and trigonometric KZ equations. 

Remark. The recent paper [CFW] studies generalized hypergeometric func- 
tions on the torus - integrals over twisted cycles of products of powers of elliptic 
functions. These functions satisfy certain linear differential equations with ellip- 
tic coefficients, but it is not clear how these equations are related to the elliptic 
r-matrix equations. 

In Section 1 we introduce a realization of an affine Lie algebra twisted by an inner 
automorphism and twisted versions of Verma modules and evaluation modules. 
This twisting is necessary to eventually produce solutions of the elliptic r-matrix 
equations. 

In Section 2 we define twisted intertwiners and deduce a differential equation 
for them. Then we define twisted correlation functions and show that they satisfy 
a twisted trigonometric KZ system. This system, however, can be reduced to the 
usual KZ system by a simple transformation. 

In Section 3 we study trace expressions of the form (6) and prove that they 
satisfy an elliptic r-matrix system of differential equations in the variables log Zj . 
(We call this system the elliptic KZ equations) Wc also deduce one more differential 
equation which expresses the derivative ^ in terms of F. Thus we get a consistent 
system of n + 1 differential equations - the extended elliptic KZ system. 

In Section 4 we show that the elliptic KZ equations are modular invariant - they 
are preserved under the action of the congruence subgroup T(N) of the modular 
group. The n + l-th equation involving the derivative by q is "almost" invariant - it 
undergoes a very minor modification under the action of an element of T(7V). This 
implies that the fundamental solution of the extended elliptic KZ system changes 
under the action of the modular group according to a certain representation of this 
group . In other words, the functions F(zi, z n \q) yield nontrivial examples of 
vector-valued automorphic (modular) forms. 

In Section 5 wc consider a few simple special cases of the extended elliptic KZ 
system for sZ 2 in which it can be solved in quadratures. In this case, solutions 
express in terms of elliptic and modular functions. Example 4 treats a slightly 
more complicated case in which the elliptic KZ system reduces to a special (but 
non-integrable) case of the so-called Heun's equation - the general second order 
linear differential equation on CP 1 with four regular singular points. 

In section 6 we compute the monodromy of the elliptic KZ equations. According 
to the results of [Chi], this monodromy yields a representation of the generalized 
braid group of the torus. We compute this representation and show that local 
monodromies (around the loci Zi = Zj) are the same as for the usual KZ equations 
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and can be described in terms of the quantum i?-matrix - a known result from the 
theory of r- matrix systems [Chi] . Global monodromies (around the T-cycle of the 
elliptic curve) are described in terms of ordered products of i?-matrices, of the type 
occuring in the quantum KZ equation (see [FR]). In the special case when for all 
j Vj is the TV-dimensional vector representation of si n ~ the monodromy represen- 
tation is a representation of the double affine Hecke algebra recently introduced by 
Cherednik[Ch3]. As an aside, the examination of monodromy helps to prove that 
if k = where M is an integer then the elliptic KZ equations are integrable in 
elliptic functions. 

A generalization of the results of this paper to the case of a quantum affine 
algebra will be described in a forthcoming paper. 
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1. A twisted realization of affine Lie algebras. 

Let g be a finite dimensional simple Lie algebra over C of rank r. Denote by 
<, > the standard invariant form on g with respect to which the longest root has 
length \pl. 

Let f) denote a Cartan subalgebra of g. The form <, > defines a natural identi- 
fication f)* — ► h: A /ia for A e f)'. We will use the notation <, > for the inner 
product in both () and 1)*. 

Let A+ be the set of positive roots of g. For a g A+, let e a , f a , H a be the stan- 
dard basis of the s^-subalgebra in g associated with a: [H a , e a ] — 2e a , [H a , f a ] — 
—2f a , [e Q , f a ] = H a . Let \a\ denote the number of summands in the decomposition 
of a in the sum of simple positive roots. 

Let N be the dual Coxeter number of g. Let p = | X) Q eA+ a - The elements p 
and hp satisfy the relations p{h a ) = a(h p ) = \a\ (note that in general H a ^ h a ). 

27ri( hp 

Let 7 be an inner automorphism of g: 7(a) = AdC(a), a g g, where C — e~~N , 
I G Z, 1 < I < N, and Z, N are coprime 1 . This automorphism is of order N. 

The action of 7 on root vectors is as follows: 7(e a ) = e'"'e a , 7(/ a ) = e _ ' a '/ Q , 
where e — e 2 ^ ll / N is a primitive AT-th root of unity. 

Let Xj, 1 < j < r, be an orthonormal basis of f) with respect to the standard 
invariant form. 

Let g = g ® C^T 1 ] ffi Cc be the affine Lie algebra associated with g. The 
commutation relations in this algebra are 

(1.1) [a(t) + Xc,b(t) +fj,c] = [a(t),b(t)] + — I < a' \t)b{t) > t' 1 dt ■ c 

27ri V |*|=i 

for any two g-valued Laurent polynomials a(t), b(t), and complex numbers A, p 
[Ka]. The elements e a ®t m , f a ®t m , Xi®t m , c, for m g Z, a g A+, form a basis 
of g. 

1 Note that throughout the paper we denote the complex number i = by a roman "i", to 

distinguish it from the subscript i, which is italic. 
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Define the subalgebra g 7 of g consisting of all expressions a(t) + Ac with the 
property a(et) = 7(a(i)). 

Lemma 1.1. (see [PS], p. .36) The Lie algebras g 7 and g are isomorphic. 

Proof. The elements e a ®tW +mN , f a <g> f-M +mJV , Xi®t mN ', c, for m £ Z, a £ A+, 
form a basis of g 7 . Define a map <p ■ 07 ~ * by: 

^>(/a®r |Q|+mW )=/a®i m , 

^•(g)f roJV ) = a; j ®t m ,m^0, 

1 c 
(1-2) 0(a;i) = ori - jjP{xi)c, 4>{c) = — . 

It is easy to check that is an isomorphism of Lie algebras. 

Thus, the twisting of g by 7 does not give us a new Lie algebra. However, the 
twisted realization g 7 of the affinization of g will be very convenient in further 
considerations. 

Let us translate some well known results about representations of g into the 
"twisted" language. 

First of all, define the polarization of § 7 : § 7 = g+ g~ © f) © Cc. Here g+ is 
the set of polynomials a(t) vanishing at 0, and g 7 is the set of polynomials a(t) 
vanishing at infinity. 

Next, define Verma modules over g 7 . This is done exactly in the same way as 
for the untwisted affine algebra. Let A £ f)* be a weight, and let k be a complex 
number. Define to be a one dimensional module over g+ © () © Cc spanned by 
a vector v such that g 7 annihilates v, and cw = kv, hv — X(h)v, h £ f). Define the 
Verma module 

Now define evaluation representations. Let V be a highest weight module over 
0. It will be convenient to assume that V is a quotient of a Verma module. Define 
the operator C £ End(V) by the conditions: Caw = j(a)Cw for any w £ V, a £ g, 
and Cwq — wq, where wq is the highest weight vector of V. 

Let V(z) denote the space of V- valued Laurent polynomials in z, and let Vc(z) 
be the space of those polynomials which satisfy the equivariance condition w{ez) — 
Cw(z). 

The natural (pointwise) action of g on V(z) restricts to an action of g 7 on Vc(z). 
For this twisted action we have an analogue of Lemma 1.1. 

Lemma 1.2. The isomorphism (f) transforms the g^-module Vc{z) into a g-module, 
isomorphic to V(z). 

Let us introduce the twisted version of currents. Set 



J ea ( z ) = e a ® t^+ mN ■ z -\°\-mN-\ 



mG2 



J fa (z) = V f a ®t-\ a \ +mN -z\ a \- mN -\ 



(1.4) J h (z) = h®t mN ■ z~ mN -\ h£ f>. 
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Thus by linearity we have defined J a (z) for any a E g. 
Define the polarization of currents: 

J+ (z) = e Q ® t^ +mN ■ z -W\~mN-i^ 

J+(Z) = fa®t-W +mN -zW- mN -\ 

m<0 

(1.5) j+(z) = ^h®i- z - 1 + Y h ® tmN -z- mN -\ ftef). 

m<0 

This defines Ja( z ) for all a6g. Now set 

(1.6) J-(z) = J+(z)-J a (z). 

Note that this polarization is not quite the same as the standard polarization of 
currents for the untwisted fl[Ka], i.e. the isomorphism <p does not match up these 
two polarizations. 

The Lie algebra g 7 has a natural Z-grading by powers of t. Thus every Verma 
module M\.k is natuarally Z-graded: the degree of the highest weight vector v 
is zero, and for every homogeneous vector w dcg(a <g> t~ m w) = deg(w) — m. Let 
d € End(M^ fc) be the grading operator: if w is a homogeneous vector then dw = 
deg(w)w. The operator d satisfies the commutation relations [d,a(t)} = ta'(t), 
[d, c] = 0. 

Let us find an expression for d in terms of elements of g - a twisted version of 
the Suguwara construction. We assume that k ^ — 1. 

Proposition 1.3. 

(1.7) 



d 



mJV 



w/iere ;: is ifte standard normal ordering: 



: e a <g> t n / Q <g> r" := 



e a ®t n f a ®t~ n n<0 
/ Q ® i""e Q <g> t" n > 



, h(E>t n h®t- n n<0 
1.8 : /i ® t n /i <g> := ^ " , /i e f). 

' h®t- n h®t n n>0 



< A, A > 
'2(fe+l)' 



Proof. Let A^a.k be the Verma module over g with highest weight A and central 
charge K. Lemma 1.1 implies that the isomorphism cp transforms the module ^ 
over g 7 into the module Ma,k over g, with K = Nk, A = A + kp. Let D be the 
grading operator in A4a.k which is associated with the grading of g by powers of 
t. Then, according to Lemma 1.1, 4>{d) — ND + h p . Therefore, 



(1.9) 



d = ^{ND + h p ) = N(j)~ 1 {D) +h p +<p,p>k 
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For the standard affine algebra g, the operator D is given by the Suguwara 
formula: 
(1.10) 

D = -kTnY, ( E : e a f»f a ®t-™:+\±: Xj ® t™ Xj t~™ = ) + ^ + ^ > , 
mGZ \aeA+ i=i / v y 

where the normal ordering is defined by (1.8) and 
(1-11) : e a f a := ^(e«/ Q + / a e a ). 

Substituting (1.10) into (1.9), after some algebra we get (1.7). 

Let us now extend the Lie algebra g 7 by adding a new element d satisfying the 
relations [d, a(t)] = ta'(t), [d, c] = 0. Denote the obtained Lie algebra by g 7 . Then 
the action of g 7 in every Verma module M\ t k extends to an action of § 7 : one sets 
d = d — 2(fc'+i) (t° S et r id °f tnc f rcc term in (1.7)). The action of d can also 
be defined in evaluation representations Vc(z): d — z-^. Thus Vc{z) becomes a 
§ 7 -module. 

2. Twisted intertwiners and Knizhnik-Zamolodchikov equations 

We will be interested in g 7 intertwining operators $(z) : M\, k — > M ly ^(g)z A Vc'(z), 
where the highest weight of V is /i, ® denotes the completed tensor product, and 
A is a complex number. 

Lemma 1.1 and the results of the untwisted theory imply the following statement: 

Proposition 2.1. Operators <f> are in one-to-one correspondence with vectors in 
V of weight X — v. This correspondence is defined by the action o/<f> at the vacuum 
level. 

Let z be a nonzero complex number. Evaluation of the operator <J>(z) at the 
point z yields an operator &(z ) : M\. k — ► M v ^ (g> V, where M denotes the 
completion of M with respect to the grading. 

From now on the notation &(z) will mean the operator <!> evaluated at the point 
z G C*. This will give us an opportunity to consider the operator $(z) as an 
analytic function of z. This analytic function will be multivalued: $(z) = z A <f>°(z), 
where $° is a single- valued function in C* , and A = <v ' V 2(k+r)' X> ■ 

Let u belong to the restricted dual module V*. Introduce the notation 
$u(z) = (1 ® u)($(z)). $ u (z) is an operator: M x . k — > M v ,k- 

The intertwining property for $(z) can be written in the form 

(2.1) [a® *"*,$„(*)] = z m <£ au (z). 

It is convenient to write the intertwining relation in terms of currents. 
Lemma 2.2. 



1 (N , N 

[Jt(Q,*u{*)] = S N ' Z ^e a u{z), a G A+; 



C \cc\-l N-\a\ 

(2.2) [Jf a (Q,*u{z)] = \ N _ CN aeA + . 
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The identities marked with + make sense if \z\ > \(\, and those marked with — 
make sense if \z\ < |£|- 

Now we are ready to write down the twisted version of the operator Knizhnik- 
Zamolodchikov (KZ) equations. 

Theorem 2.3. The operator function $ u {z) satisfies the differential equation 

(* + 1)JU„(*) = E ( J l( Z )®f a u(z) - ®f a u(z)J-(z)) + 

ueA+ 

(2.3) 

r 

q£A+ j=l 

Proof. The logic of deduction is the same as in the untwisted case (see, e.g.,[FR]). 
Equation (2.3) is nothing else but the intertwining relation between $(z) and d: 

(2.4) z^* u (*) = -[0, $„(*)]■ 

Substituting the expression for d (Eq. (1.7)) into this relation, after some calcula- 
tions we obtain (2.4). 

Now let us define the twisted correlation functions. Let V\,...,Vn be highest 
weight representations of g, and let $ % Ui {zi) : M\ uk — * M\ i _ 1 .k, 1 < i < n be 
intertwining operators. Set A„ = A and Ao = v. Let v\ is the highest weight vector 
of M\ t k, and let v* be the lowest weight vector of the restricted dual module to 
M Vj k- Consider the correlation function 

(2.5) $> Uu ... tUn (z 1 ,...,Z n )=<vt,$ 1 Ul (z 1 )...$Z n (Zn)v\>, Ui&Vi, 

which makes sense in the region \z\\ > ■ ■ ■ > \z n \. The function ^ can be regarded 
as taking values in the space V\ (g> • • • <g> V n . 

The function if satisfies a twisted version of the trigonometric KZ equations. 
Before writing these equations down, let us recall the definition of the trigonometric 
classical r- matrix [BeDr] . 

Let fl = J2aeA+( e a®fa+fu® e °')+J2j=i Xj®Xj, £1 e 000. For < p < N-l, let 
P be the eigenspace of 7 in with the eigenvalue e 2mp / N . Let be the orthogonal 
projection of £1 to the subspace g p ® 0_ p of £g> g. The trigonometric r-matrix has 
the form 

(2-6) r(z) = f + E 

p=0 

Introduce the convenient notation: 

a t = Idy 1 <g>...<g>v i _ M ®a® lv i+1 ®-®v n , a e U(g), l<i<n; 

(2.7) ifr = Ea S ® &S J a s ,b s eU(g), then r<j = J^a")^'),-. 

The main property of the trigonometric r-matrix r(z) is the classical Yang-Baxter 
equation: 

(2.8) [rij(zi/zj),rjk(zj/zk)] + [rjk{z /z k ),r k i{z k /zi)] + [r ki (z k / z i ),r lj (z l / zj)] = 0. 
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Theorem 2.4. The function \f satisfies the following system of differential equa- 
tions: 

(2.9) {k + l)z i -^- = Y,r iJ {z i /z J )* + -{hx + h v )V. 

1 j¥=i 



Proof. The proof is analogous to that in the untwisted case (see [FR]). It is based 
on the direct use of relation (2.3). 

Since the twisted intertwining operators are obtained from the usual ones by a 
simple transformation, we should expect that system (2.9) should reduce to the 
untwisted trigonometric KZ equations. This turns out to be the case. Indeed, set 
Q = zf, and *(&,.. .,Cn) = {z" lp )i...{z l n p ) n ^{z ll ...,z n ). Then the function 
satisfies the trigonometric KZ equations: 

(2-10) N(k + l) Ci § = £ n ^ + "^ g + hh x + h» + 2h p )*, 

where are the half Casimir operators defined in the introduction, and A = X+kp, 
v = v + kp. Therefore, the standard theory of the KZ equations can be applied to 
the study of the properties of (2.9). 

3. Traces of intertwiners and elliptic r-matrices 

From now on the letter g will denote the Lie algebra sIn(C) of traceless N x N 
matrices with complex entries. The dual Coxeter number of this algebra is N, and 
the rank is N — 1. The Cartan subalgebra [) is the subalgebra of diagonal matrices, 
and the element C is the matrix diag(l, e~ 2 , ...,s~ N+1 ) (up to a factor). 

Let B be the N x N matrix of zeros and ones corresponding to the cyclic per- 
mutation (12. ..TV). Note that BC = eCB. 

Define a new inner automorphism [3 of g: /3(a) = BaB -1 , a e g. This automor- 
phism has order TV and commutes with 7: (3 o 7 = 7 o (3. 

The action of the automorphism B naturally extends to g and g 7 . Note that on 
g 7 , unlike g and g, B is an outer automorphism: it corresponds to the rotation of 
the affine Dynkin diagram of g (which is a regular iV-gon) through the angle 2-k/N. 

The action of (3 in g 7 preserves degree, hence, it preserves the polarization. 
Therefore, it transforms Verma modules into Verma modules. In other words, we 
can regard B as an operaator B : ~ > M^x),ki where by convention [3(\)(h) = 
A(/3 _1 (ft,)). This operator intertwines the usual action of g 7 and the action twisted 
by [3: [3(a)Bw = Baw, a € g 7 , w e M\^. 

Let v — (3~ 1 (X), and let & 3 u .( z j) be as above (cf. section 2). Let q be a complex 
number, < \q\ < 1. Following the idea of Frenkel, Reshctikhin ([FR], Remark 2.3) 
and Bernard [Bcr], introduce the following function: 

(3.1) F Uli ... iUn ( Zl , z n \q) = Tr \ Mv:k «(*!) . . . <S>ZJzn)Bq- d ). 

This function takes values in Vi<8> • • • ®V/v. From now on it will be the main object 
of our study. 
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It turns out that the n-point trace F(zi, ...,z n \q) defined by (3.1) satisfies a re- 
markable system of differential equations involving elliptic solutions of the classical 
Yang-Baxter equation for sIn- Let us deduce these equations. The idea of the 
method of deduction is due to Frenkel and Rcshctikhin. 

Differentiating by Zj , we get 

d 

(k + l)—F(z u ...,z n \q) = 

Tr ($1 (zi) ■ ■ • (Zj) ~ & faUj (zj)J~ (zj)) • • • K n [zn)Bq- d ^ + 

Tr ^< (zi) ■ ■ ■ Jt (zj)H aUj (Zj) H aUj (z 3 )J u ( Zj )) . . . 3> n Un (z^Bq-^ + 
(3.2) 

Tr (*IM) ' ' ' E W^Xu^Zj) - n iUl (zi)J^(zj)) ■ ■ ■ KJ^Bq-^ . 

Now let us pull the currents J + , J~ around: the currents J~ will move to the 
right up to the end, then jump at the beginning and continue to move to the right, 
and so on; the currents J + will move to the left up to the beginning, then jump at 
the end and continue to move to the left, and so on. As we do these permutations, 
we will need relations (2.2) and also the following identities: 

(3.3) q - d jt(z)=qjt(qz)q- d , 

(3.4) BJt(z) = J± a) {z)B. 

After J + and J~ have made M complete circles, equation (3.2) will have the 
form 

d 

(k + l)—F(z u ...,z n \q) = 
TrUi i (z 1 )---^(Q M j;M {e Jq M z j )^ aUj (z j )- 

®luM)<l- M Je-M {e Jq- M z J )) ■ ■ ■ KJza)Bq- a y 
Tr(<(zi) • • • ]>> M JpM {f Jq M z 3 m aUj (zj)- 

& eaUj (z 3 )q- M Jp- M{f Jq- M z 3 )) ■ ■ ■ KSz n )Bq- d ^ + 
, N-l 

TrK>) ■ ■ ■ ^ {q M J+ M{xi) {z j )<^ XiUj {z j )- 
^ l=i 

H l u j (z j ^~ Mj 0-^^ q ~ Mz ^ ■ ■ ■ K n (z n )Bq- d y 

n 

(3.5) ^XgFizu...,^), 



i=i 
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where 

(3.6) 

and 



EM yiv, 
p=-M ij 



M.p 



spm x M,p . _ . 

l^p=-M,p^O ^jj 1 ~ J 



^ = £ 



,p A N-l-\a\M 



qP N Z? - zf 



-(P p (e a Mf a h 



(3.7) 



£ 



qP(qP Zi ) 



\a\-l 



qP N z N _ Z N 



-(^(/ Q )) 4 (e Q ),+E 



1 (qP Zi ) N + z N 



2 Zl qPN z N _ Z N 



Now we want to pass to the limit M — > oo. Right now we cannot do so since the 
limit does not exist. In order to be able to pass to the limit, we should write down 
equation (3.5) for M = L, L + 1, L + N— 1, add these N equations together, and 
divide by N: 

d 

(k + l)—F( Zl ,..., Zn \q) = 

L + N-lr , 

jr E Tr <(^)---E^ M 4«(e a) (^) $ i« 3 fe)- 

M=L L ^ a 
$ i« 3 fe)9" M ^"-M(e a )(9" M ^)) • • • KJZa)Bq- d y 

<« 3 .fe)9- M ^--M (/a) (?- M ^)) • • • *Z n {zn)Bq-°y 

, N-l 

^ 1=1 
KuM^ M Jp-M {xi) (q- M z,)) ■ ■ ■ KJzn)Bq- d ) + 



(3i 



£* 

1=1 



F(2i,...,z n |g), 



In the obtained equation it will already be possible to take the limit. Moreover, 
since J2p=o P P (h) = for ft € fj, the part of the right hand side of (3.8) involving 
currents associated with the Cartan subalgebra elements will disappear as M — > oo. 
The same will happen to the currents associated with the root elements e a , f a 
because these currents do not contain a term of degree —1 which is the only term 
that could possibly have given a nonzero limit. Thus, in the limit we get a simple 
equation: 



(3.9) 



_d_ 

d Zq 



F(z u z n \q) = E X « 
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where X99 = lim^ % EmJI" 1 X% . 

The function X°° admits a simple description in terms of elliptic functions. 
First of all, it is easy to check that X^ = for any M. 

Next, looking at the function Xf- can be represented in the form Xfj = j-pff [z%l Zj), 
where p M ' (z) is a rational function with values in q ® q: 



M 



p \N-\a\ 



p M (*)= E 



p=-M VgA + 
„, ■ , ,| N - 1 



(3.10) Yl Q pN Z N 1 

Therefore, X?? = —p%j{zi/zj\q), where p(z\q) is an elliptic function of log z with 
values in g (g) g (we have used notation (2.7)). We can tell what this function looks 
like by looking at its residues. 

From (3.10) we see that the only poles of p(z\q) are at the points e m q p , and all 
these poles are simple. The residue of p(z\q) at z = e m q p is equal to that of p M {z) 
for M > \p\, i.e. it equals 



Res \ z = emqP p(z\q) = E (^ mH /a®/ ?P (0 + 



(3.11) e m ^e a ® (F(f a ) + £ x ; ® /F(z*)) = ^ (1 ® 7 " m /3 p )(«) 

Because of the obvious homogeneity property J ^2 i Zi^- = 0, we have the unitarity 
relation Pij(z) = —pji{—z). 
Let q = e 27riT . Let 



(3.12) ((x\r) = -+ lim ]T 



1 



+ y- 



x- j - It (j + It) 



-M<j,l<M,p+P>0 

be the standard Weierstrass function. Equation (3.11) implies that 



(3.13) 
1 

27riiV 



E 



(i® 7 - ro w) 



/ log z m 

~ N ~ 



P t\Nt ) - C( -jj -pr\T 



0<m,p<N -l,m 2 +n 2 >0 

Thus, we have proved 

Theorem 3.1. The function F(z\, z n \q) satisfies the system of differential equa- 
tions 



(3.14) 



dF 



(k + l) Zi — = J2p ij (z i /z j \q)F, 1 < i < 



n. 
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The function p(z\q) is a solution of the classical Yang-Baxter equation. It is easy 
to see that this function is nothing else but the elliptic r- matrix for sIn due to 
A.Belavin [Be]. Thus, we have given a representation-theoretical interpretation of 
the local system associated with the elliptic solutions of the classical Yang-Baxter 
equation. 

Remark. A.Belavin and V.Drinfcld [BcDr] showed that elliptic solution exist 
only for the simple Lie algebra sIn, and every nondegenerate elliptic solution is 
equivalent (=conjugatc) to const • p(z\q) (for a suitable primitive N-th root of unity 
s). 

Observe that equations (3.14) transform into the KZ equations (2.9) as q — > 0. 

_ <A,A> 

This was to be expected since lim 9 ^ <7 2 < fc+1 > F(z\q) = ^f(z). We will call equations 
(3.14) the elliptic Knizhnik-Zamolodchikov (KZ) equations. 

There is one more equation satisfied by the function F - a differential equation 
involving the first derivative by q. Differentiating F by q, we obtain 

(3-15) - q — = TT \ Mxk (& Ui ( Zl )...<f> Un ( ZN )Bq- d d). 

Plugging the expression for d in (3.15), we see that in order to obtain the differential 
equation for F, it would suffice to express the traces 

Tr \u Kk « (zi) • . . <S> n Un (z N )Bq- 9 f a ® t-^- mN e a ® t^+ mN ), 
Tr \ Mx , k « (zi) . . . $Z n (z N )Bq- d e a ® t^~ mN f a ® 
(3.16) Tr \ Mx k « ( Zl ) . . . (z N )Bq- 9 xi ® ^ m % ® t mN ) 

in terms of the original trace F(z\q). This can be done as follows. 

Take the expression for the first trace in (3.16) and move the factor e Q ® t' Q ' +mW 
from left to right. When it has made N full circles, we will have 

Tr \ Mx , k « (zi) . . . <S>Z n (z N )Bq- d f a ® t-H- mAr e Q ® = 

n N-l 

J2 z \f+rnN q . p{lal+mN)Tr w ^ $i_ p(ea) ^ (z .) . . . ( ZAr)i?g -9 /Q t -|a|- m i V) + 

j=l p=0 

q -N(\ al+mN ) Tl ($ i i ^ (^)Bg- 8 (/i a + (M + mA )fc))+ 

(3.17) 

It is easy to see that for any h £ f) 

(3.18) Tr | M ,, fc « (*i) (^)B«- S ft) = -(1 - r 1 )" 1 W^^k). 

Similarly we have 
(3.19) 

n 

Tr Uf A , fc «(z 1 )...^„(z Ar )i? g - 9 .f Q ^-l«'- mAr ) = -^^ |a| - mJV (l-/3- 1 3 |a|+roJV )- 1 (/a)^(z|?). 
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Thus, we obtain the following expression for the first trace in (3.16): 

Tr \ Mx , k « (zi) . . . K n (z N )Bq' a f a ® r^~ mN e a ® *M+™ W ) = 

n 

^ ^lal-m^al+mJV^ _ ^M+^-l (/q) . (1 _ ^-M-r^-l ^ jF{Aq) + 

(3.20) (1 - g -"(M+»»"))-i ^ ( i _ ^- 1 )" 1 (-/ la ) + (|a| + mAT)^F(z| g ). 

The second and the third trace in (3.16) are treated quite similarly, and finally 
after some calculations we obtain: 



(3.21) 



where 



t,3 — l J 



q)F, 



L(z\q) = Yl 2~ H ~ miV (l - /rV a|+TnJV ) _1 (/a) ® (1 - P^q-^-^yHeah 

aeA+ m>0 

^2 J2 z |a| - roJV (i - /3- i g- |Q|+mW )- i (e a ) ® (i - f3- 1 q lal ~ mN y 1 (.f a )+ 

a£A+ m>0 
JV-1 

^ J2 z- mN (l - Pr x q mN Y\x p ) ® (1 - /3- 1 g- mA, )- 1 (x p ) + 

p=l m>0 

2V-1 

(3.22) - J2 a - r 1 )- 1 ^) ® (i - r 1 )- 1 ^) 

is a function with values in g ® g (once again, we use notation (2.7)). Note that 
Lij = Lji. Hi — j, then by La we mean p,(L)i, where fi : U(q) <8> U(q) — ► C/(g) is 
multiplication: /u(a ® 6) = a6. 

The form of equation (3.21) can be simplified. Indeed, the consistency of (3.21) 
and (3.14) for all k implies that 

q ^=z±(Liz\ q ) + Hz-*\ q) ), i.e. 



(3.23) L(z\q) + L(z- 1 \q) = 2L(l,q) + J 



qdp(z\q) 



dz. 



z dq 

Therefore, we finally get the theorem (in the formulation we use notation 2.7)): 
Theorem 3.2. The function F satisfies the differential equation 



dF 



q) F, 

ij 



(3.24) 



.<*>=/! 



qd p(z\q) 
8q 



dz 
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Corollary. 

(3.25) Tr \ MQtk (Bq- d ) = 1. 

Proof. Apply Theorem 3.2 to the case of a single intertwining operator &(z) : 
Mo,k — * M 0: k ® Vc( z )' wnere ^° is the trivial representation of g. In this case it is 
obvious that $(z) = Id. Therefore, (3.23) simply follows from the fact that (3.22) 
acts by zero in V°. 



4. Modular invariance of the elliptic KZ equations 

Let us now discuss the modular invariance of the elliptic KZ equations. Introduce 
new variables yj — — ^fp 1 , t = From now on let us use the notation k = k + 1. 
Rewrite equations (3.14) and (3.24) in the new coordinates: 

(4-1) k— = Y,P*M - Vj\r) F > 1 < * < n 

(4.2) - E 



where 



dr 

1,3 = 1 ]<i 



P*(y\r) = ^p(z\q), 



s*(y\r) = 2^ S {z\q) = J" ^^dx, 

(4.3) L*(t) = 27riL(l|g). 

Let T(N) be the congruence subgroup in 5*L 2 (Z) consisting of the matrices equal 
to the identity modulo N. We have the following almost obvious property. 

Proposition 4.1. Equations (4-1) are invariant with respect to the group T(N). 
That is, they are preserved under the change of variables jji = f = if 

Proof. Partial derivatives with respect to the new coordinates are given by 

_d_ _ (cT + d) _d_ 

djji dyi ' 

(4.4) |: = (- + ^ + + %*|t- 

Thus we only have to prove that p*(y\f) — (cr + d) p* (y\r) . To do this, it is enough 
to observe that both sides of this equation 

(i) have simple poles at the points of the lattice m + pr with the same residues 
i(l®7- m /3 p )(«), and 

(ii) satisfy the unitarity condition: p\ 2 {y\ T ) — ~/°2i( — v\ T )- 



17 



Now let us study the invariance properties of equation (4.2). First of all, it is 
easy to see that 

(4.5) = (cr + d f 9 J^l + c(cr + rff^i. 
ot or ay 

Integrating this equation against dy, we get 

(4.6) s*(y,f) = (cT + d) 2 s*(y\T)+c(cT + d)(^yp%y\T)-Yimxp*(x\T)y 

By our definition, lim^o xp*(x\r) — ^. Thus, under the change of variable 
t — > f, y — > y equation (4.2) transforms into the following equation: 

n2 ^ nV- 
(CT + d) K— + C(CT + d)2_^ ^-0^: = 

(4.7) 

]T(ct + d)V( yi - w |r)F + ]T c(cr + d) f(y t - yj )p*( yi - y \r) - ^\ F + £ L«(f)F. 

Combining (4.7) with (4.1), we can get rid of the derivatives by yt and reduce (4.7) 
to the form: 

( 4 - 8 ) 4^ = E U fob-) - t^w) ^ + («■ + d r 2 E 

Now we need to find the law of transformation of L*(t). 
Lemma. 



(4.9) L*(f) - (cT + d) A L*(T) + 



2N 



Proof. Let C(r) = L*(f) - (cr + d) 2 L*{r) - c{cT 2 + N d)n e <8> 0. We know that both 
equations (4.2) and (4.8) are consistent with (4.1). Therefore, we have 

(4.10) [Cn(r) + 2C ia (r) + C 22 (r) + ^ + ^ + ^ , P* 12 (y\r)] = 0. 

Observe that fin = A ® 1, f2 22 = 1 ® A, where As [7(g) is the Casimir element. 
Since the Casimir element commutes with the Lie algebra action, (4.10) reduces to 
the relation 

(4.11) [Cn(T) + 2Ci2(r) + C 2 2(T),pI 2 (y|r)]=0. 

This relation has to hold for all y, which implies that the expression C(t) = Cn (r)+ 
2Ci 2 (t) + C 22 (t) commutes with all 7-invariant elements in g (g) g, i.e with all 
elements of the form Xi <g> Xj, e a £g> (3 p (f a ) and f a <g> /3 p (e a ). 
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We are going to prove that C(r) = 0. Let 

(4.12) C(t) - < T ) Xi ® ^ + S M T X e « ® Wa) + W«) ® e a ). 

Pick arbitrary two elements X and Y in rj. Consider the expression 

c xy = [x ® i, [i ® y, 2C]] = [x ® i, [i ® y, <5]] = 

(4.13) -2 £ & ap (a(X)a(/TP (F))e Q ® /?"(/«) + a (Y) a (/r"pO)/ Q ® /3 p (e Q ) 

(for brevity we do not specify explicitly the dependence on r) . This expression has 
to commute with Xi ® Xj for all This immediately implies that Cxy = for 
all X, Y, i.e. b ap = for all a, p. Therefore, 



(4.14) C = J2 a *- x * 



> Xi. 



Let Xi = 1 <g> Xi + Xi <g> 1. Then C* = a i^ 2 - Therefore, 

[C^® /?"(/„)] = 
(4.15) Yl - ^)( I i))«i( I i e « ® W«) + e « ® Wa)*0- 

i 

In order for (4.15) to be zero for any a,p, we must have ^ a ((l~ f3~ p )(xi))a,iXi = 
for all a, p, which can only happen when all a, are zero (because the roots span 
the space f)*) Q.E.D. 

The Lemma we just proved implies 

Proposition 4.2. Under the change of variables (y, r) — ► (y, f) system of equations 
(4-1), (4-2) transforms into a system of equations equivalent to the combination of 
(4-1) and the following equation: 

i,j=l j<i i—1 

where Aj is i/ie Casimir operator in the i—th factor of the tensor product ■ ■ ■ ®V n . 

Thus, the system of equations (4.1), (4.2) is almost invariant under T(N): the 
first n equations are unchanged under the action of this group whereas the last 
equation gets a very simple extra term. 

Consider the fundamental solution T{y\r) of the system (4.1), (4.2) - a solution 
with values in End(Vi® . . . ®V n ) determined by the initial condition 

(4.18) T(y\, ...,y n \r) — > Id as r — > ioo, Imj/j+i — Imy., — > +00 for 1 < j < n — 1. 
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Theorem 4.3. (on the modular invariance of solutions of the elliptic KZ equa- 



(4.19) r{y u ...,y n \f) = (cr + d)^ A \F(t/!, ...,y„|r)x(A), 

where x(A) is an operator in ■ ■ - ®V n dependent only on the 2x2 matrix 

A (i.e. independent of yi andr). 

Assume that the representations Vj are irreducible. Then Aj are simply complex 
numbers. In this case the function x(A) is a projective representation of T(N): 

(4.20) X(A 1 ) X (A 2 ) = a(A 1 ,A 2 ) X (A 1 A 2 ). 

The 2-cocycle a(A\,A2) is very easy to describe. Let A = \ ^ ) € 5X 2 (Z) 



x c d / 

be called positive if cither c > or if c = but d > 0; otherwise, let A be called 
negative. Clearly A is positive if and only if —A is negative. Let \A\ denote A if A 
is positive and —A if A is negative. Then 

{1, 1^4.! 1 1 A 2 1 is positive 

e 2« K ; |Ai||y4 2 |is negative 

Of course, this cocycle is a coboundary since H 2 (T(N),Z) = 0, so the action of 
T(N) in the projectivization of Vi(g>V^(g> . . . ®V n comes from a linear action. 

Thus, the theory of the elliptic KZ equations gives us a natural method to assign 
to every set of irreducible finite dimensional representations V\,...,V n of sIn an 
action of the congruence subgroup T(N) of SX 2 (Z) in the tensor product of these 
representations, V\ <g> • • • <g> V n . 

In fact this construction allows us to obtain a representation of the entire mod- 
ular group SL 2 {Z). For this we need to assume that Vi are finite dimensional 
representations of GLjv(C) for all i. For brevity we will also assume that N is odd 
(this assumption is not very essential, but in the even case one has to be a little bit 
more careful). 

We will need to use the Weil representation of the group SL^Z/NZ). This rep- 
resentation is defined as follows. Take the iV-dimensional vector space U — C N and 
define an action of the Heisenberg group H N =< x, y\x N = y N = 1, xyx~ 1 y~ 1 commutes with x, y > 
in this space by x — > B, y — > C. This is the basic irreducible representation of Hn- 

Now observe that SL 2 (Z/NZ) acts by automorphisms of H^: if A = ( a ^ ] then 



A(x) = x a y b , A(y) = x c y d . Moreover, the representation U A of Hm obtained from 
U by twisting of U by A is isomorphic to U . Therefore, the group SL 2 (Z/ NZ) pro- 
jectively acts in U in such a way that Azu = A(z)Au, A e SL2(Z/NZ), z e Hn, 
u e PU 7 and this action is unique. The space U with the constructed projective 
action of SL 2 (Z/NZ) is referred to as the Weil representation. This representa- 
tion defines a homomorphism W : SL 2 (Z/NZ) — > PGL N (C). Since the group 
SL 2 (Z/NZ) for odd N does not have nontrivial central extensions, this homomor- 
phism lifts to a map W : SL 2 (Z/NZ) -> GL N (C). 
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Proposition 4.4. Let A = d) ^ SL 2 (Z), and let the change of variables 

(y,r) — ^ (y,T~) be defined as in Proposition 1^.1. Then 

(4.22) f(y u .., y n \f) = (cr + d)^ ^^>0(A)F(y u ...,y n \r) X (A), 

where x(A) is a projective representation of SL 2 (Z) in V\ <g> • • • <g> Vjv, and 9 is the 
composition of the three maps: 

Aut{Vi <g> • • • <g> V N ), 



(4.23) 5i 2 (Z) SL 2 (Z/NZ) GL N (C) 



where itj : GLjv — » ^.Mi(Vj) are i/ie homomorphisms defining the action of the group 
GL N in Vi. 

The proof of this statement is simple and similar to the proof of Theorem 4.3. 
The Weil representation of SL 2 (Z/NZ) arizes naturally when we consider the trans- 
formations of the lattice L of poles of p*(y\r) which do not preserve the lattice of 
periods NL. 

The function x(A) satisfies equation (4.20) with the 2-cocycle a still being defined 
by (4.21) (now for the entire SL 2 (Z)). Since H 2 (SL 2 (Z), Z) = 0, this cocycle is a 
coboundary, so x(A) comes from a linear action of SL 2 (Z). 

Remark. It is not clear how to compute the representation X (A) for any 
nontrivial example. A good example to start with would be N = 2, n = 1, and V\ 
is the 4-dimensional irreducible representation of sl 2 . In this case it seems that V 
will be a direct sum of two irreducible 2-dimcnsional representations of T(2), and 
the solutions of (3.24) will be some nontrivial vector-valued modular functions. 

5. Some examples 



Let us consider the special case N = 2, g = sl 2 . In this case (3 acts as follows: 



/3(e) - /, (3(f) 
(5.1) 
where 



e, 0(h) 



-h. Therefore, we have 



p( z W) = 2 



a(z\q)(e <g> e + / <g> /) + b(z\q)(e <g> / + / <g> e) + c(z\q)h <g> h 



a(z\q) = — 



logz 
27ri 



logz 
27ri 



2r 



logz 1 
2th 
log z 



~ 2 -r\2r 



2iri 



t\2t 



+ ao(q), 

Mi), 



(5.2) 

c(z\q) 



1 

47ri 



logz 
27ri 1 



2r +C 



log 2 

2th 



r|2r 



flog. 



c 



log z 
27ri 



x - r|2r 



and the constants a , bo, c are chosen to satisfy the condition p(z\q) = — p(— z\q). 
We also get from (3.22) 



£(l|?) = -E 



o 2m+1 (l + 9 4m+2 ) 



m>0 



(1 - g 4m + 2 ) 2 



2g 



Am+2 



(5.3) 



m>0 
1 /l 

2 



~4m+212 



(ec 



ie + /« 
)e)+ 



/) 



^ (1 + g 2m ) 2 

m>0 v y 7 



h®h. 



co(q), 
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Now let us calculate some nontrivial 1-point traces. 

Example 1. Let $(z) : M_\^ — * M\.k <8> Vq{z) be an intertwining operator, 
where V 1 is the two-dimensional irreducible representation of 5/2- Then we must 
have A = ±5. Denote the corresponding operators by $ ± , and introduce the 
notation: T±(q) = Tr($ ± B(7 _a ). It is clear that these traces do not depend on z. 
Let us compute them. 

We have proved that T± (q) satisfy the equation 

(5.4) K q^=L(l\q)F. 

In the case we are considering the matrix L turns out to be a scalar 2x2 matrix: 

(5.5) L(l,q) | vl = + £ m)2 -E (1 _ 4m+ 2 )2 - 

m>0 v y 7 m>0 V * ' 

Therefore, we can explicitly integrate equation (5.4), which gives us the following 
answer. 

Let v± be the basis of V , such that hv± = ±i>±, ev_ = v+, fv + = V-, ev + = 
fv- = 0. Let T)(q) be the Dedekind function: 

00 

(5.6) V (q) = q 1 ^ J] (! 1 m )- 
Then 

(5.7) T±(q) = V (q 2 )£v±. 

Since the Dedekind function is a modular form of weight 1/2, the function (5.7) 
is a modular function of weight 3/8/t, which is by no means a surprise in view of 
formula (4.16). 

Example 2. Let &(z) : M-\ t h — > Ma,*: <8> V^c(z) be an intertwining operator, 
where l^ 2 is the three-dimensional irreducible representation of s/2. Then we must 
have A = or ±1. Denote the corresponding operators by $°, $ ± , and introduce 
the notation: T (q) = Tr(^°Bq~ 9 ), T±(q) = Tr($ ± Sg- a ). These traces are com- 
puted similarly to Example 1. The matrix L(l, q) is now no longer scalar, but it is 
a diagonal matrix, so equation (5.4) is still easy to integrate. Here is the answer: 

Let w + ,«o,w- be the basis of V 2 , such that hv + = 2v +1 hv- = — 2u_, hvo = 0, 
ev_ = fv + = vq. Then: 



(5.8)T ± ( g ) = ^( g 4 )^(^ 



T (q)=r 1 (q A )-r 1 (q 4 )--v Ql 



Let us now calculate a simplest 2-point trace (a part of this computation is due 
to A.Kirillov Jr.). 
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Example 3. Consider intertwining operators: <& (z) : M\,k — * M^ii,* ® ^cK- 2 )- 
We can combine four traces out of these operators: 

(5.9) T ±± (z\q) = Tr(* ± (z 1 )$ ± (z 2 )Bq- d ),z= -. 

It is possible to calculate these traces explicitly using the elliptic KZ equations. 
We have proved that the traces (5.9) satisfy the equation 

dF 

(5.10) k— =p(z\q)F 

In the case under consideration, the traces take values in the four dimensional 
space Vq ® Vq, and this equation can be explicitly solved. Indeed, let us seek the 
solution in the form 
(5.11) 

g( z \q) = g++(z\q)v+®v + +g + -(z\q)v + ®v-+g- + (z\q)v-®v + + g — {z\q))v- 
Consider the functions 



hi(z\q) = 


\(g++(z\q) + g- 


-(z,q)), 


h 2 (z\q) = 


\{g++{z\q) -9- 


-(«»«)), 


h 3 (z\q) = 


\{g+-{z\q) + g- 


-+{z,q)), 


h 4 (z\q) = 


\{g+-{z\q) - g- 





(5.12) 

We have 
(5.13) 

g(z\q) = h 1 (z\q)(v ++ -V—)+h 2 (z\q){v ++ +v — )+h 3 (z\q)(v + _-v_ + )+h 4 (z\q)(v + _+v. 

Equation (5.10) yields a separate first order linear differential equation for each 
of the functions hi, hi, I13, h 4 : 

K ^z- = 2 { - a + c)kl > 

K ^z- = 2 {a + C)k2 > 
«-^ = 5 (-6-c)ft3, 

(5.14) K ^ = I( 6 _ c ) fe4 . 

These equations are easily solved. To write down the solutions, it is convenient 
to use the function 

(5.15) ~p'{v\2t)= Yl (y-m-2pry 
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In terms of this function, the solutions of (5.14) can be written in the form: 

h 1 (z\q) = C 1 {q)E(-z\q)&, 
h 2 (z\q) = C 2 (q)E(-qz\q)^ , 

h 3 (z\q) = C 3 (q)E(z\q)^ , 
hi(z\q) = C 4 (q)E(qz\q)^ , 

(5.16) 

The coefficients Cj(q) are easily found from equation (3.24): 

c 1 {q) = c 1 r } ^)Her-(^)\ 

c 2 (q) = cM^v(<i 2 r-(^)\ 
CM = c 3 , 

C±{q) = C i r ] {q 2 t^tf)-i, 

(5.17) 

It remains to say what values of constants Cj correspond to the traces (5.9). 
This information is given below: 

T ++ : C\ = C 2 = 1, C3 = C4 = 0; 
T — : & = -1, C 2 = 1, C 3 = C 4 = 0; 
T_+ : Ci = C 2 = 0, C 3 = -1, C 4 = 1; 

(5.18) T+_ : d = C 2 = 0, C 3 = C 4 = 1. 

In general, solutions of the elliptic KZ equations cannot be expressed in terms of 
classical elliptic and modular functions. Their components are more complicated 
special functions associated with an elliptic curve. 

Example 4. Consider the elliptic KZ equations with coefficients in Vq <S> V£- 
Let i>i, V-i be the basis of Vq introduced in Example 1 (earlier we used the notation 
v+,V- for this basis), and let W- 3 , W-i, w\, w 3 be a basis of Vq in which fw 3 = 
wi, fwi = 2w-i, fw-i = 3w_ 3 , fw- 3 = 0, ew- 3 = io_i, ew_i = 2wi, ewi = 
3u>3, ew 3 — . Let us look for solutions of the elliptic KZ of the form: 

F(z\q) = h + (z\q)(vi ® Wi + V-\ ® W-\) + h-(z\q)(v-i ® W-i — V\ ® 

(5.19) f+(z\q)(vi <g> w_3 + u-i <S> W3) + ® w_ 3 - t)_i <g> 103). 

The functions /i± and /± can be found from the following 2x2 linear systems of 
differential equations: 

K ^7 = ( 2 C + a) ^ ++ 2 5/+ ' 

(5.20) ^ = ^ + - ^c/ + ; 
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K ^T = { 2 C - a)h - + 2 bf -> 

(5.21) K ^ = hh_-Lf_. 
y ' dz 2 2 

Let us reduce the systems (5.20) and (5.21) to second order differential equations. 
It will be convenient to use the notation <f>± = ± a. Differentiating the first 
equation in (5.20) and (5.21), we get 

(5.22) Kh'l - <j>' ± h± + 4>±ti ± + hf ± + h'f ± . 
Substituting the second equation in (5.20), (5.21) into (5.22), we obtain 

(5.23) nh'l = 4>' ± h ± + cj>±ti ± + i^b(^bh± - |c/±) + \b'f±. 

Substituting the relation f± = f — 4>h±) into (5.23), and taking into account 
the fact that —b'/b = 2c, we deduce the second order equation for h: 

(5.24) h'i + ^(2 + Oc T « _1 a) ft± - ^ (<t>'± + |;& 2 + (|: + 2 ) c< ^±) ft ± = °- 

Equation (5.24) is a second order equation on the elliptic curve C/ < l,2r > 
with four regular singularities at the points of order 2. This equation is obviously 
invariant under the transformation z — ► —z. Therefore, if we make a change of 
variable w — p(z\2r) (here p{z\2r) = — Q'{z\2t) is the Weierstrass p-function), 
we will obtain a second order equation on CP 1 with rational coefficients and four 
regular singularities, at the points E\ = p(l/2|2r), E 2 = p(r|2r), E 3 — p(r + 
1/2 1 2t), and oo. By a linear transformation of the independent variable the singular 
points can be placed at 0, 1, A, and oo, where A = ^ 2 Z%\ ■ 

A general second order equation with four regular singularities was first studied 
by K.Hcun[Hc] in 1889, and it is now called Heun's equation. The general form of 
Heun's equation is 

^ ^ ' (.t x — 1 x — \) x(x — l)(x — A) 

where 

(5.26) a + /J-7-(5-e+l = 0. 

Its solutions (which are called Heun's functions) were studied in much detail in the 
middle of the 20-th century [Er]. 

A Heun's equation is characterized by a P-symbol which contains the information 
about the eigenvalues of monodromy matrices at singular points (see [Er]). The P- 
symbol is a table whose columns correspond to the singular points of the equation. 
In each column there are three numbers - the coordinate of the singular point and 
the two roots of the characteristic equation at this point. The P-symbol of (5.25) 
is 



{0 1 A oo n 

a z \. 
1-7 1-5 l-e (3 > 
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Of course, the P-symbol does not determine the equation uniquely since it contains 
no information about the constant p. 

In the special case of equation (5.24), we can easily write down the P-symbol. 
The P-symbol of the transformed equation (5.25) for the plus sign is 

{0 1 A oo n 

-£ + § z\. 
Ill L 4. 3 J 

2 k k 2k ' 4 

For the minus sign the columns under and A in (5.29) have to be interchanged. 

It is easy to see that equations (5.20) and (5.21) (and hence (5.24)) cannot 
be solved in elementary and elliptic functions. To verify this fact it is enough 
to look at the limit q — > O.In this limit we are supposed to get the solution of 
the KZ equations in the space V 1 <£> V 3 . This solution is known to be given by a 
hypergeometric function. The combination of parameters in this specific case is such 
that this hypergeometric function does not express through elementary functions. 
This implies that solutions of systems (5.20), (5.21) cannot be expressed in terms 
of elliptic functions. 

6. Monodromy of the elliptic KZ equations 

In this section we will study the monodromy of the elliptic KZ equations with re- 
spect to the lattice of periods, and compute the monodromy matrices. Although it 
is difficult to compute the solutions, the calculation of monodromy is fairly straight- 
forward. 

Let us first describe how to interchange the order of intertwining operators. 

Let $ w > A ' I/ ( / j) : M\,k — * M Vj k(&z A Vc(z) be the intertwining operator such that 
< v*, Q w ' ' v (z)v\ >= w, w € V . Suppose that 2:1,22 are nonzero complex 
numbers, and we have a product $ w i< x i> x °( Zl )$ W2 ' X2 > Xl (z 2 ) ■ M\ 2 ,k —> Mx . k <g> 
V\ (g> V 2 where V\ and V 2 are finite dimensional representations of g. The question 
is: can this product be expressed in terms of products of the form ^{z2)^{z\)l Of 
course, we can only talk about such an expression after analytic continuation, since 
the former is defined for \z\\ > \z2\, and the latter for \z\\ < \z2\- However, if we 
apply analytic continuation, the answer to the question is positive, and given by 
the following theorem. 

Theorem 6.1. (see [FR],[TK]j Let be a basis ofVi~ x °, and let be a basis 

f V ^2-U Then 

(6.1) 

$xrx 1 MM( Zl )q>v.* 1 ,)*M( Z2 ) = A ± Y J R%rsX l u{ x 2AQ) VlV2 ^ Vi "' v ^°{z2)^ XivM ' v {z 1 ), 

where is the analytic continuation along a path in which z\ passes 22 from the 
right (for plus) and from the left (for minus), respectively, P ± (A, [j,) VlV2 is a matrix, 
and a is the permutation: V\ ® V2 — ► V2 ® Vi ■ 

Clearly, the matrix P ± (A, n) VlV2 represents a linear operator 
(Vi O V 2 -» (V 2 ® Vi) x -^. Therefore, if we define 

(6.2) P±(A) yiV2 = @ tl R ± (\ii) VlV2 , 



(5.29) 
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then R ± (X) VlV2 will correspond to an operator: V\ ® V2 — > Vi <S> V\ . This operator 
has been computed [Koh] , [Dr] , [TK] , [SV] , and it turned out to be proportional to the 
product of the quantum i?-matrix of the quantum group Uq(g) and the permutation 
a (the order of this product depends on the sign, plus or minus), where q — e 2m / NK . 
Hereafter we will assume that the matrix .R ± (A) ViV2 is known. 

Let us now compute the monodromy of solutions of the elliptic KZ equations. 
Again, we will need to assume that V\, V n are finite dimensional representations 
of GL N . 

The elliptic KZ system is a local system with singularities, £, on the space 
E n , where E is the elliptic curve C*/T, and V is the multiplicative subgroup in 
C* generated by q N . The fiber of this local system is V\ ® • • • ® V n . Such an 
interpretation, however, is not very convenient for computation of monodromy since 
C has too many singularities: they occur whenever Z{j ' Zj = e m q p . It would be more 
natural to regard the elliptic KZ system as a local system on the n-th power of a 
smaller elliptic curve E — C*/T, where T is generated by q and e. In this case, the 
singularities would occur only on the loci Zi — zj. But unfortunately, the elliptic 
KZ system is not a local system with singularities on E n : its right hand side is not 
q, e-periodic. Therefore we would like to produce a local system on E n starting with 
C. For this purpose we will use the fact that C has a finite group of symmetries. 

Recall the notation of Section 4: Hn denotes the Heisenberg group of order N 3 . 
Let H N denote the n-th Cartesian power of H^. Since the group is naturally 
embedded in GLn (it is generated by B and C), we have a natural representation 
of Hm in Vi and hence a representation of H' N in V\ ® • • • ® V n . On the other hand, 
the group Hm operates on E: Bz = qz, Cz = e~ 1 z, so the group naturally 
operates on E n . These two actions can be combined into an action of H^j in the 
trivial bundle over Em with fiber V\ ® • • • ® V n . This action has the following 
property. 

Lemma. The group H N preserves the local system C 

Proof. The lemma follows from the definition of the elliptic r-matrix (formula 
(3.13)). 

Now we can create a new "local system" S = C/H N . The fiber of this local 
system is no longer the space V\ ® ■ ■ ■ ®V n but rather the quotient of this space 
by the action of the center of H N (which is, of course, not a vector space). In 
this section we will describe the monodromy of this local system. This monodromy 
will be a linear representation of a suitable central extension of the fundamental 
group of i?™\{diagonals}, and it obviously contains all the information about the 
monodromy of the elliptic KZ equations. 

As before, we will use the fundamental solution T{z\, z n \q) defined by (4.18). 
This solution takes values in the space End(Vi<g)- • -<g)14.). The fundamental solution 
has a defining property: if u is a vector in V\ ® • • • ® V n then Tu is the solution of 
the system (3.14), (3.24) which tends to u as q — > and Zi/ Zi + \ — > 00. 

Below we study the monodromy of the fundamental solution. This is equivalent 
to studying the monodromy of the local system S. 

We will use the representation of E as a parallelogram on the complex plane: 
= jf + yT, < x, y < 1, and write z = {x,y). We choose a base point 
(z?, z°) = « ...,x° n ,y°) on E n such that 1 > x\ > ■ ■ ■ > x° n > 0, 1 > y? > 

■ • • > y° n > 0. 
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It follows from the definition of T{z\, z n \q) that it can be represented in the 
form z^z® 2 ■ ■ ■ z^ n J r o{z\, z n \q) where D\,.. ., D n are operators in Vi ® • • • <g> V n 
defined as follows. Let vectors Vj € Vj satisfy the condition hvj = Xj(h)vj, h £ f), 
Xj € f)*. Let 



(6-3) A, = (/3-l)- 1 (^ x ,) + E^' 0<j<n. 

i=i i=i 

Then 

(6.4) Dj(v! <g> • • • <g> u n ) = - Aj,Aj > ~ < Aj_1 ' Aj_1 > ■ «i <g> • • • <g> u„. 

(from this definition of Dj it immediately follows that [Di, Dj] = for any 

This observation helps us find the monodromy of the function T around the 
e-cycles on E n . If Zj is rotated around the origin anticlockwise, through the angle 
2iril/N, with the rest of the variables fixed, the function T multiplies by the matrix 
e 2-KiWj/N f rom the right, and undergoes a conjugation by C^-the action of C in Vj. 
We denote these monodromy operators by Ej : 

(6.4) Ej{F) = CT x TCje^ iW *l N . 

The monodromy of the fundamental solution in the neighborhood of the locus 
Zj = Zj+i can be found with the help of Theorem 6.1. Using this theorem and the 
representation of solutions as traces of products of intertwiners, we immediately 
find that the monodromy of the function J 7 , around the locus Zj — Zj+± (in the 
anticlockwise direction) is 

(6.5) hi+iiF) = F s j,j+i, 
where 

S jJ+1 =Rj(V u ..., V j+1 , Vj,..., V^R+iVu . . . Vj,V j+1 , V n ), 
and the linear operators 

Rf{Vi, V n ) : Vt ® • • • ® Vj ® Vj+x ® • • • ® V n -» Vi ® • • • ® V j+1 ® Vj ® • • • ® V n 
are defined as follows: if Vi e Vi, 1 < i < n, and hvi = Xi(h)vi, Xi € i)* , h G f), then 
fi^(Vi, K)(«i«>- • ■0^®Wj+i®- • •«>«„) = ui<8>- • ^^(Aj+i)^^ 1 ^®^!)®- • -C 

where Aj are defined by (6.3). 

This monodromy is obviously the same as for the usual (trigonometric) KZ 
equations, whose solutions are given by matrix elements of intertwiners rather than 
traces (see section 2) - a known fact which was first proved by I.Cherednik[Ch2]. 

Let us now find the monodromy of the function T around the g-cycles. Consider 
the cycle on E n in which Zj passes Zj + \, z n from the right, hits the circle \z\ = \q\ 
(which is identified with the circle \z\ = 1 through the map z qz), jumps over to 
the circle \z\ — 1, and then passes z\, ...,Zj-\ from the left, returning to its initial 
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position (this corresponds to the path yj (t) = y® — t mod 1 , with the rest of Xi and 
yi remaining unchanged). Using Theorem 6.1 and the expression of T in terms of 
traces, we find the monodromy matrices Qj for the described g-cycles. Indeed, we 
have to interchange &(zj) in the trace expression with $(zj + i),...,$(z„), then with 
q~ d and B, and then with $(zi),...,$(zj_i). This results in an expression of Qj 
as a product of i?-matrices. To write this expression down, introduce the following 
notation: if s e S n is a permutation of n items then set (s) — (V s ^ , . . . , V s ( n ) . 
Let tj be the elementary transpositions + 1), and let Sj m = t m -\ . . .tj + \tj, 
j < m < n, Sj m = t m . . . tj-2tj-i, 1 < m < j (we make a convention that for two 
permutations <J\,<J2 ^l^ti) — ci(c2(i)), 1 < j < n, i.e. the factors in a product of 
permutations are applied from right to left). Then the operators Qj are expressed 
as follows: 
(6.6) 

Qi(F) = B 3 TR+(lA)R+ +1 ( 3 , J]+l ) . . . fl+_ 1 ( Sj „_ 1 )B7 1 Er( Sjl ) . . . i^fe-i), 

where Bj denotes the action of B in Vj. 

Remarks. 1. Expressions similar to (6.6) (ordered products of i?-matrices) 
occur in the theory of correlation functions for quantum affine algebras. Such corre- 
lation functions satisfy a quantum analogue of the Knizhnik-Zamolodchikov equa- 
tions - a system of difference equations discovered by Frenkel and Reshetikhin[FR] . 
The structure of the right hand side of this system resembles (6.6). It is not clear 
if it is merely a coincidence or not. 

2. It is seen from the definitions of Ej and Qj that [Ei,Ej] = [Qi,Qj] = for 
any 

As we have already remarked, the monodromy of the local system S defines a 
representation of a central extension of the fundamental group of the complement 
of the diagonals Zi — Zj in E n - the pure braid group of the torus. To describe this 
representation in more detail, let us assume that Vi = V for all i, where V is some 
finite-dimensional GLjv-module. This does not cause any loss of generality since 
we can always set V = V\ © • • • © V n to include the previously considered case. But 
now the elliptic KZ system (and the local system S) has an additional symmetry 
- the symmetry under the simultaneous permutation of the variables zi and the 
factors in the product V <g> V ® • • • ® V. Therefore, the monodromy representation 
can in fact be regarded as a representation of a certain (in general, not central) 
extension of the full braid group of the torus. 

The braid group of the torus, BT n , is generated by the elements T i} 1 < i < n— 1, 
Xi,Yi, satisfying the defining relations 

TiT i+1 Tt = T i+1 TiT i+1 ; T{Tj = TjT it j > i + 1; 
(TiXi) 2 = (X^f: (T^- 1 ) 2 = (l^Ti) 2 ; 

TjX x = XtTj, TjY\ - Y\Tj, j > 1; 
X 2 Y- 1 X- 1 Y 1 - If, where X l+1 = T^T,; 
(6.7) X Fi = YtXo, where X = X X X 2 ...X N 

It is also convenient to define the elements Yj + \ = T^ l YjT~ l and Y = Y{Y 2 . . . Y„ 
To picture the braid group of the torus geometrically, one should imagine n 
"beetles" crawling on the surface of the torus C/(Z + rZ) starting from some fixed 
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positions zi,...,z n (zj = N~ l Xj + ryj, x\ > x 2 > ■■■ > x n , yi > y 2 > ••■ > y n ) 
so that at no time two beetles can be at the same point and after some period 
of time (say 1) the "beetles" return to their original positions, possibly with some 
permutation. Then the beetles will trace out some collection of curves in Torus x 
[0, 1] - a braid diagram. Such diagrams can be composed by attaching the bottom 
of one of them to the top of the other. Under this composition law, braid diagrams 
form a group - the braid group of the torus BT n . The element Tj corresponds to 
the intertwining of the j-th and j + 1-th braids (the j-th and j + 1-th "beetles" 
switch, the j-th "beetle" passing the j + 1-th one from the right), and the elements 
Xj and Yj arise when the j-th "beetle" crawls around the x-cycle and y-cycle of 
the torus, respectively, in the negative direction of the x-axis (respectively, y-axis), 
with the rest of the "beetles" unmoved. 

Now we are in a position to formulate the result about the monodromy of the 
local system S. 

Theorem 6.2. The monodromy representation of an extension of BT n associated 
to the local system S is defined as follows: 

X, i — ► e-^^Cr 1 , 

Y i ^ ■ ■ ■ Rn-\ B j ■ ■ ■ RJ-i 

(6.8) Tj » (R+y 1 . 

Remark. Since all the spaces Vj are the same, we have dropped the permuta- 
tions labeling the i?-matrices. 

Proof. The theorem follows from Theorem 6.1 and formulas (6.4)-(6.6). 

The extension of BT n involved in Theorem 6.2 is very easy to describe. Consider 
the quotient of the group BT n by the relations Tf = 1. Then we will obtain the 
group S n K (Z™ © Z"), where S n acts on both copies of Z™ by permutations of 
components. Let H be the Heisenberg group - the central extension of Z © Z by 
Z by means of the 2-cocycle w(x, y) = xiy 2 — x 2 yi- Then we can construct the 
group S n k H n (H n is the n-th Cartesian power of H) which is a rank n abelian 
extension of S n x (Z™ © Z"). Let us denote the pullback of this abelian extension 
to BT n by BT n . The monodromy representation of S is thus a representation of 
BT n in V ® • • • O V. 

Now assume that V is irreducible. Then instead of the group BT n we will 
have to deal with a rank 1 central extension BT n of BT n which is constructed as 
follows. Take the group S n K (Z"©Z n ), construct a rank 1 central extension of this 
group by adjoining a new central element c satisfying the relations XiYi = YiXiC, 
XiYj = YjXi, i 7^ j, and then pull this extension back to BT n . It follows then that 
the monodromy representation of S is a projective representation of BT n which 
comes from a linear representation of BT n . 

The group BT n is closely related to the double afhne braid group of type A„_i 
defined by I.Cherednik in his recent paper [Ch3] - the group generated by the 
elements Ti,X i} Yi and a central element 5 satisfying modified relations (6.7): the 
relation X 2 Y^X^ 1 Y 1 = T? has to be replaced by X 2 Y^ 1 X^ 1 Y 1 = 5T?. Let us 
describe the connection between them. 

Observe that the elliptic KZ system commutes with the diagonal action of the 
Heisenberg group -ff^ lag generated by X n and Y in V <£> • • • <£> V. Decompose 
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V ® • • • ® V into a sum of irreducible representations of H N m& : V ® • • • ® V — 
®j Pi <2)Wi, where Pi are distinct irreducible representations of Hn, and Wi are 
multiplicity spaces. Then each summand Pi <g> VFi is a subrepresentation of -BT„. 
Let <fo : -BT„ — > End(Pi <g> Wj) be the corresponding homomorphism. 

Assume that n and iV are coprime. Then we can define a homomorphism £ : 
ST„ -» H n by = C 1 /", = S 1 /™, £(T;) = 1. (Here 1/n is regarded as 

an element of the ring Z/7VZ.) Composing this homomorphism with the diagonal 
action of Hn in V (g> ■ ■ ■ (g> V, we get a projective representation V'iO) of 5T„ 
in End(Pi <g> Wi). Notice that ipi(g) — fa(g) if g € H N ias . Let us write </>j as 
01 (5) = V'i(s) ' ^i(g)^ 1 4'i(g)- Because the elliptic KZ system commutes with H^ ag , 
ipi(g) commutes with ip i (g)~ 1 fa(g) , which implies that ipi(g)~ 1 fa(g) = Id ® Xi(sO> 
where x«(') is some projective action of BT n in Wj. We also have ipi{g) = V'iCfl')® 1> 
where ipi is an action of the group in Pi. Therefore, we have fa = ipi®Xi- Therefore, 
we can easily compute the 2-cocycle on BT n corresponding to \i as the difference 
of the 2-cocycles for fa and tp. This cocycle is the pullback from S n k (Z™ © Z") of 
the 2-cocycle given by: 

n 1 n 

^((^.xSy 1 ), ( S2 ,x 2 ,y 2 )) = YMlvf - xU) - 

1 ^ ■ 1 

Z— 1 J — 1 

(6.9) 3i,2 €5„, x^y'^Z". 

The extension of BT n by means of this cocycle is exactly the double affine braid 
group. We denote this group by BT n . 

The special case when V is the TV-dimensional vector representation of GLn is 
especially interesting. In this case, we are getting a representation of the double 
affine Hecke algebra This algebra was recently defined by I.Chcrednik [Ch3] 
as the quotient of the group algebra C[BT n ] by the relations (Tj — q)(Tj + q^ 1 ) = 
0. Indeed, the matrix is diagonalizable and has the eigenvalues q 1 ^^ and 

— q _1_ 2iv j (recall q = e 2nl / NK ). This statement easily follows from the fact that 
this matrix is a monodromy matrix of the elliptic KZ equations. Therefore, the 
matrix Rj = q^R^ satisfies the equation 

(6.10) (R+-q)(R+ + q~ 1 ) = 0. 
Thus, the correspondence 

(6.11) X 1 1 ^ e-^ UDl / N C^\ Yi ^ R+R+ . . . fl+^Bf 1 , Tj ^ R+ 

defines a projective representation of BT n in the space V <g> V ® • • • <g> V which 
can be written as P® W, where P is the TV-dimensional irreducible projective 
representation of BT n obtained by composing the homomorphism £ defined above 
with the standard action of Hn in C^, and W is a representation of the double 
affine Hecke algebra Sj^. 

The element S acts in the representation W by multiplication by 
Remark. In fact, the monodromy of the elliptic KZ system (3.14) extended 
by equation (3.24) yields a representation of the semidirect product SZ^Z) x BT n 
in the space End(F <g> • • • <£> V). As we have already remarked in section 4, it is 
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not clear how to compute the modular part of this monodromy, i.e. the action of 
SL 2 (Z) . This computation, at least for one nontrivial example, is a very interesting 
and challenging problem. 

As a conclusion, let us note that the study of monodromy helps us to find out for 
what special values of k the elliptic KZ equations are integrable in elliptic functions. 

Proposition 6.3. If k = l/MN, where M is an integer, then the matrix elements 
of the fundamental solution T(z\, z n \q) are finite products of rational powers of 
theta functions of expessions q m e p Zi/ ' Zj, < m, n < N — 1. 



Idea of Proof. If n = 1/NM then q 
the permutation of factors, so 
(s is rational). Then we have 



= 1. Therefore, R + is a scalar matrix times 
= ax Id, and bjj+i{T) = aT . Let a = e 27T[s 



(6.14) F(z u ...,z n \q)=l[ [] e(q m ePz l /z ] \q N ye^( / uj 

i<j 0<m,p<N-l ^ 

where w is a matrix- valued elliptic differential form on E n and 

(6.15) Q(z\q) = Y[ (1 q m z)(l - q^z- 1 )^ - q m+1 ). 

m>0 



Remark. A similar result holds for the conventional KZ equations (1): if 
k + h v = 1/M, where M is an integer then solutions of the KZ equations are 
algebraic functions. 
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